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We demonstrate a technique based on noise measurements which can be utilized to 
study dynamical processes in protein assembly. Direct visualization of dynamics in 
membrane protein system such as bacteriorhodopsin (bR) upon photostimulation are 
quite challenging. bR represents a model system where the stimulus-triggered struc- 
tural dynamics and biological functions are directly correlated. Our method utilizes 
a pump-probe near field microscopy method in the transmission mode and involves 
analyzing the transmittance fluctuations from a finite size of molecular assembly. 
Probability density distributions indicating the effects of finite size and statistical 
correlations appear as a characteristic frequency distribution in the noise spectra 
of bR whose origin can be traced to photocycle kinetics. Valuable insight into the 
molecular processes were obtained from the noise studies of bR and its mutant D96N 
as a function of external parameters such as temperature, humidity or presence of an 
additional pump source. 
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I. INTRODUCTION 

Noise features from a functioning system can reveal valuable insight into processes oc- 
curring over several time scales^. Stochastic fluctuations are observed in systems like gene 
expression, protein, nano-medicine and can be used to probe and characterize gene circuits 
and phase transitions^. Identifying sources of noise, controlling and processing noise-related 
data using sophisticated techniques based on frequency domain analysis provides insight 
into molecular information that control the system^"-. Noise as a tool to study electronic 
events and processes at microscopic length scales is extensively utilized in condensed matter 
physics. There are not many situations encountered where noise studies of optical absorp- 
tion processes at room temperature has proven to be useful. Noise in biological processes 
are dealt at macroscopic-systems level and their correlation to molecular events are rarely 
emphasized. It is expected that fluctuations in a small system to be pronounced and we uti- 
lize that aspect by choosing a model biomolecular organization with the appropriate system 
parameters. Typical fluctuations from finite sized and correlated systems yield probability 
density functions (PDF) that deviate from Gaussian or Lorentzian distributions^. Kinetics 
associated with biochemical reactions can be modeled based on modified Langevin equations 
(LEQ). Solution of LEQ expressed in terms of reaction rates for mRNA and protein systems 
has been used to obtain noise frequency ranges which in turn are related to the decay rates 9 . 
It has also been shown that fluctuations in small systems are larger in magnitude and can 
get comparable to the mean values of the variables^. Hence there is a need to modify the rate 
equations to accurately model finite systems. Effective mesoscopic rate equations (EMRE) 
have been formulated to account for the breakdown in the linear noise approximation^"—. 
Bacteriorhodopsin (bR), a retinal protein is an ideal system to demonstrate the utility of 
noise studies, as it can be modeled by a 2-state model LEQ in spite of a complex photo- 
cycle. The longest time constants are in the ms-s range and can be controlled by external 
parameters. Finite or small system sizes are essentially defined by the film thickness and the 
near-field scanning optical microscope (NSOM) probe geometry controlling the excitation 
volume. 

NSOM based sub-diffraction limit imaging technique utilizes narrow optical fiber tips to 
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image single molecules and nano-systems^. It was shown in our laboratory, that it was 
possible to monitor the intermediate states of bR, using transmission NSOM 14 . The light 
driven proton pump mechanism in bR is initiated by the isomerization of retinal protein 
from all-trans to 13-cis configuration followed by formation of series of intermediate states, 
referred to as J, K, L, M, N and O states^ 1 ^. M state (absorption maxima at wavelength, 
A ~ 412 nm) representing the deprotonated state has the highest population among all 
other intermediate states and exhibits strongest spectral shift following photo-excitation at 
A ~ 570 nm which finally decays thermally or upon blue light illumination. The distinct 
spectroscopic signature of the intermediate states in the photocycle of bR provides the 
platform to observe fluctuations associated with the different molecular states. 

To our knowledge, the intensity fluctuations time series in the transmittance of light 
through a molecular assembly has not been analyzed and related to internal molecular 
changes. Although single molecule spectroscopic techniques have been demonstrated for 
fluorescent systems^. A probable reason could be from the experimental constraint of 
non-availability of systems with appropriate inherent signal to noise ratio within the data 
acquisition rates. The advent of NSOM techniques overcomes some of these constraints, 
and in the process signals from an assembly consisting of < 600 molecules can be closely 
studied with a combination of probe and pump beams (SN.l)^. A distinct time series 
signal overriding on a noise feature of sizable amplitude signal describes the transmission 
fluctuations Tr(t). The characteristic signature of bR in form of pump induced absorption 
is observed to accompany changes in Tr(t). Tr(t) upon transformation to the frequency 
domain provides a consistent picture of the events. We implement an algorithm to process 
the data acquired over large time scales to arrive at a robust representation of the frequency 
distribution. 

II. THEORY 

The kinetics involved in the reaction initiated by the photoexcitaion of the ground state of 
bR is reasonably well understood. The photocycle kinetics, primarily the M-state lifetimes in 
bR are known to be affected by pH conditions, additional light-pump at A405 corresponding 
to the M-state excitation, humidity, temperature and presence of metal nanoparticles^. 
These changes introduced by external factors are manifested as peak and linewidth shifts 
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in the noise frequency profile. A quantitative handle on the problem can be initiated using 
a simplified model where the bR molecule is reduced to a system with two excited states 
(Fig. [Tti), where B, M and N represents the ground state, excited state, and a long-lived 
intermediate respectively with a finite reversible rate along with the associated rate constants 
k{ (B — > M), k\(M — > B) for the i th molecule. k{ = a(X)I, where a(X) is the absorption 
coefficient and I is intensity. k\ on the other hand is the inverse of the M state lifetime 
(tm)- A set of LEQ can be formulated for the two state model with noise factor as an 
additional term to justify the nondeterministic outcome: dNs/dt = —k{ Nb + + <t>B{t) 

and dN M /dt = k{ N B — k'-N M + (j) M (t) where, N B and N M are concentrations of B and M 
state of a similar set of i molecules and 4>b and 4>m are noise terms. The nature of terms 4>b 
and 4>m could be white noise with or without correlations. It should be mentioned that the 
choice of Ornstein Uhlenbeck process with correlations can lead to more realistic solutions. 
In the presence of A 40 5 for exciting the M-intermediate state, these equations get further 
modified associated with the additional pump-rate. The noise terms in these equations 
can have a characteristic noise frequency range analogous to the approach by Simpson et 
al^. The PDF can now be attributed to predominantly arise from M state lifetime and can 
assume a sum of Lorentzian profiles (modified by normal distribution). 

In finite systems the trimers in bR are known not to get simultaneously excited, due to 
heterogeneity— (Fig. [Tb). The fluctuations arising from these sources get further modified 
by the cooperative effects, random excitation and decay of molecules, reversible states in the 
photocycle and thermally driven transitions. Heterogeneity is averaged out in bulk systems, 
but in a finite size ensemble they play a significant role in extracting useful molecular 
information^. The general principles in the method may be applicable to other analogous 
chromophore containing membrane-protein systems. 

III. MATERIALS AND METHODS 

Wild-type bR (WT-bR) films of different thickness were obtained from aqueous suspen- 
sion of purple membrane patches (~ 0.5 mg/ml, pH ~ 9.2)(SF.1,2)^>^. These films were 
characterized by absorption spectroscopy and transient absorption to ascertain bR features 
(SF.3) 18 . After the films were prepared on quartz substrates the topography (AFM) (SF.4)^ 
and corresponding transmission (NSOM) map of the sample segment with a probe laser (532 
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nm) were simultaneously recorded (SF.5)^. Once the bR patch was located, the tip was 
then positioned onto a point on the patch using software controlled stepper motors and the 
piezo stage. Typically bR patches have lateral dimensions ranging from 500 nm to 2 fim 
(for multilayer) and heights ranging from 7 nm (monolayer) to several 200 nm (30 layers). 
The NSOM tip (100 nm diameter) illuminates only a fraction of the region of the 500 nm 
to 2 fim bR patch. The photoexcitaion volume is governed by the aeral coverage of the 
NSOM tip and the thickness of the patch. A typical Tr(t) trace for WT-bR (Fig. [2^ (inset)) 
was obtained using a customized Nanonics MultiView 4000 NSOM model with cantilever 
tip geometry (SF.6) 18 . In a typical experiment, a continuous wave laser (probe) was cou- 
pled to the NSOM tip (~ 100 nm) and the transmitted optical signal from a bR patch was 
measured using a photo-multiplier tube (PMT) through the objective with magnification 
50x and N.A = 0.45 (SF.7)^. Additional pump beam (405 nm) could be introduced in the 
set-up and the measurements with appropriate filters were also carried out (SF.8,9)^. Noise 
histograms for bare quartz also do not reveal any distinct features when additional pump 
source is introduced (SF.10)^. The noise histograms for the quartz substrate are primarily 
a measure of source and detector fluctuations. 

Temperature variation was carried out by connecting conducting adhesive-tapes from the 
scanning stage to a hot plate. Increasing the hot plate temperature resulted in heating the 
sample stage in the temperature range 296 K to 318 K. The temperature at the sample was 
calibrated with a digital thermometer. A scan was carried out at room temperature and 
the transmission data was recorded. Subsequent scans were carried out as the temperature 
was increased. Humidity variation was carried out by enclosing the scanning stage in a 
chamber with a steam inlet (SF.ll) 18 . A source of steam was connected and the scan was 
then carried out. The humidity of the enclosure was independently verified and controlled 
by the duration of the steam flow into the chamber. Scans were typically carried out in 
steps of 30 minutes. 

Algorithms used to analyze data were similar to the ones used in photon counting ex- 
periments to determine the coherence times of pseudo-thermal sources (SN.2) 23 . Data from 
PMT was captured using a digital oscilloscope (LeCroy, Wave Runner 6100A) at 250 kHz 
sampling rate for WT-bR and 500 Hz for D96N mutant. About 100 consecutive data sets of 
1 s (for WT-bR) and 10 s (for D96N) were captured as separate windows. For the WT-bR 
case each window was further split in 15 sub-windows and Fast Fourier Transform (FFT) 
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was applied to each one of these windows. The frequency corresponding to the maximum 
amplitude in the FFT signal was stored for all the windows. A histogram was constructed 
of the frequencies and their occurrences leaving out the DC contribution. 

IV. RESULTS 

These traces are utilized to obtain noise histograms using a rigorous, standardized proce- 
dure. Autocorrelation function (ACF) obtained from Tr(t) for quartz and WT-bR samples 
of different patch heights, clearly indicate the absence of correlations for the bare quartz 
and a sizable correlation for the bR region (Fig. Wfr)- The frequency distribution arrived 
from Tr(t) for WT-bR patch was consistently in the form of a lognormal type distribution 
about a characteristic 0J max . Lognormal distribution refers to a skewed distribution where 
higher frequencies extend the contribution, and whose natural logarithm follows a normal 
distribution given by the expression, f(x;/i,a) = (1/xaJ (2n))(exp(—(ln(x) — /i) 2 /(2a 2 )), 
where /i and a are the mean and standard deviation of ln(x) respectively. A common feature 
arrived upon Tr(t) analysis is the presence of characteristic maximum in the amplitude cor- 
responding to a frequency, uj max riding a distribution (Fig. [2b, c). The magnitude of 0J max ~ 
500 - 700 Hz is in the range of the lifetime associated with the M-state which constitutes a 
large fraction of the entire molecular photocycle span. The interpretation of the dominant 
frequency associated with the intensity fluctuations to the bR photocycle is consistent with 
a set of measurements where photocycle rates were intentionally affected by varying certain 
external factors. The overall profile of the noise spectrum is then expected largely to be 
controlled by a distribution arising from the heterogeneity in the molecular state in the 
ensemble and the photocycle kinetics. 

The skewed nature or the non-normal probability density function (PDF) can also point 
to the statistically correlated nature of these fluctuations (Fig. |2b,c). Trimers in bR are 
known to effects arising from correlations in excitation and de-excitation processes 2 ^. It has 
been shown that the degree of skewness can be related to the correlation or to the system 
size 8 . It is difficult from our results to quantify the relative contributions from these factors. 
We speculate that both correlations and finite size effects appear to play an important role 
in giving rise to the non-Gaussian PDFs which defines our system. 

We do not discount other possible representations such as combination of Lorentzians 
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and Gaussians, but a straight forward lognormal fit yielded reasonable fit parameters. The 
appropriateness of lognormal function is apparent in form of good Gaussian fits when the 
distribution is plotted on a log frequency axis (Fig. [2b, c inset). Further Kolmogorov- 
Smirnov tests were done to confirm that the statistics is lognormal (ST. 1,2)^. The lognormal 
behaviour of the distribution is an indicator of statistically correlated events or networks of 
interacting elements in the present finite size system^. It arises from multiplicative product 
of many independent random variables. A shift of u max was observed when the additional 
photo-excitation pump source (A = 405 nm) is introduced in the Tr(t) measurements. It is 
known that life time of the M state is altered by the thermally driven events in photocycle 
hastened upon the pumping the intermediate states, thereby biasing the population to the 
ground state. The blue shift in the oo max upon A 40 5 excitation can then be associated with the 
reduction in the photocycle span. The parameters in the distribution are also dependent on 
the number of bR layers (film thickness) (Fig. [2b, c). The number of bR layers is a direct way 
to control the ensemble size. The Tr(t) data exhibits a maximum at 505 Hz for 40 nm thick 
patch and 630 Hz for 100 nm thick patch which gets shifted to 630 Hz (Fig. [2b) and 720 Hz 
(Fig. [2b) upon additional photo-excitation (ST.l)^. Changes in the mean value upon fitting 
to a lognormal profile, where mean takes up the value of m = exp(fi + cr 2 /2), are similar to 
the changes in u max as a function of thickness and pump excitation. It is observed that 25 
nm to 150 nm thick patches illuminated by 100 nm tip yields PDFs which reflect internal 
molecular, M state dependence on external conditions. In case of monolayer which consists 
of about 200 trimers in the cross section of the incident light beam where quantifiable near- 
field absorption is observed, the characteristic noise pattern however is missing (SF.9)^. 
The noise features apparently appear to be dominated by the laser fluctuations and detector 
characteristics (similar to bare quartz substrates). On the other hand, in the regime of thick 
films (> 200 nm), heterogeneity and other random molecular events are averaged out and 
the noise pattern is again featureless^. Hence there seems to be a critical size regime (14 
nm - 150 nm) which yields characteristic stochastic features. 

The T dependent bR photocycle kinetics is modified in dried-film form by the presence 
of competing factor of moisture concentration related to the proton uptake process . Upon 
heating the bR films (to 320 K), the increase in moisture depletion with T reduces the 
photocycle rate. These changes are observed in the form of the red-shift in the u m with 
increasing T (Fig. [3^, inset). However, upon increasing the moisture concentration, the 
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expected increase in the photocycle rate in form of the blue shift of oo m is observed from 
lognormal distribution fits (Fig. [3b, inset). It is interesting to note that the time series data, 
Tr(t) is also a gauge for ambient conditions, as clearly noted from measurements where the 
proton concentration gets depleted and enriched upon switching off and on the source for 
moisture respectively. In order to generalize the procedure to other systems we demonstrate 
the method using genetically modified bR (variant D96N), which has a longer M state 
lifetime (> 100 ms in film) due to the replacement of aspartic acid (Asp-96) by asparagine 
(Asn-96)^. The noise measurements were suitably extended to the low frequency range and 
the data was collected at 250 KHz and 500 Hz sampling rates. Noise features in D96N: (i) 
Noise distribution around the maximum of 120 mHz was observed, which shifted to 480 mHz 
upon the additional A 40 5 excitation (Fig. (3t). A blue shift was also observed upon increasing 
the humidity (Fig. ^jp inset) (SF. 12)^. (ii) A frequency distribution was also observed about 
750 Hz which however is independent of the illumination. This can correspond to the rise 
fluctuations which are of the order of ms and dependent only on intensity of illumination 
and the absorption coefficient (SF.13)^. Low frequency (< 2 Hz) histograms for WT-bR do 
not exhibit distinct pattern or shifts with pump excitation (SF.14)^. 

The Tr(t) measurements on bR films provide an important toolkit as a model system to 
examine the role and implications of fluctuations. Signature of characteristic noise which 
appears to be ensemble size dependent in the light-induced proton pump function throws 
open interesting questions. Molecular processes involved in the functioning of the protein are 
expected to be dependent on the organization prevailing at much higher length scales and 
the collective behaviour which emerges depends on these interactions. The interactions can 
further be probed by varying external parameters and the lognormal distributions analysis 
of Tr(t) can help quantify the collective response. The general principles in the method 
should be applicable to other interesting biological systems. 
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FIG. 1. Simplified two level energy level scheme of bR showing variability introduced by different 
sources of noise, (b) Hexagonal arrangement of bR trimers depicting a possible configuration 
suggesting the statistical source for randomness 
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FIG. 2. (a) Autocorrelation functions for quartz and various bR patches, (inset) Typical time 
series noise data. Noise power spectra histograms for (b) 40 nm WT-bR patch, (c) 100 nm WT-bR 
patch, (blue bars, 532 nm illumination, red bars, 532 nm along with pump illumination). (Insets) 
Corresponding histograms plotted in log frequency scale (blue dots, 532 nm illumination, red dots, 
532 nm along with pump illumination,) and Gaussian fits (blue curve, 532 nm illumination, red 
curve, 532 nm along with pump illumination) 
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FIG. 3. (a) Red shift of noise power spectra histograms as a function of temperature for WT-bR. 
(Inset) Temperature dependence of mean noise frequency for WT-bR. (b) Blue shift of histograms 
as a function of humidity for WT-bR. (Inset) Humidity dependence of the mean noise frequency 
(H1<H2<H3<H4). (c) Noise frequency histograms for D96N mutant (green bars, 532 nm illumi- 
nation, blue bars, 532 nm along with pump illumination). (Inset) Humidity dependence of peak 
noise frequency for D96N mutant showing variation with humidity.) 
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